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The argon spectral function is constructed and applied to neutrino-argon cross section computa¬ 
tions in the plane wave impulse approximation with the Pauli blocking final state interaction effect 
taken into account. The approximations of the construction method are critically analyzed using 
the example of oxygen for which more detailed computations are available. An effective description 
of nucleus based on the information contained in a spectral function is proposed. It is demonstrated 
that its predictions are close to those obtained from the complete spectral function. The effective 
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I. INTRODUCTION 

In recent years, a lot of interest has been concen¬ 
trated on improving the description of neutrino inter¬ 
actions [ij. The motivations behind this research are 
future precise neutrino experiments which require better 
knowledge about the cross sections. It also makes crucial 
whether the theoretical models can be implemented in 
Monte Carlo (MC) event generators. 

Many papers focus on neutrino energies of a few GeV 
because this is the proposed energy range of several long- 
baseline experiments. In this energy region one has to 
consider both quasielastic and inelastic charge current 
(CC) processes. It is common to factorize neutrino- 
nucleus interaction into two steps: an interaction with 
a single bound nucleon (impulse approximation) and 
then reinteractions inside nucleus. The main difficulty 
is caused by the need of precise evaluation of the nuclear 
effects. 

Existing MC codes rely on various versions of the Fermi 
gas (EG) model because it is straightforward to imple¬ 
ment them in the MC routine. The simplest EG model 
with two adjustable parameters—Fermi momentum and 
average binding energy 0— can be modified by inclusion 
of the recoil nucleus energy in the kinematics Q and by 
local density effects It is, however, clear from the 
experience gained in describing the electron-nucleus in¬ 
teractions that the model can be considered only as the 
first approximation and that the MC simulations based 
on the EG model can be neither detailed nor precise. 

The approach which introduces more realistic picture 
of a nucleus is the one using the spectral function (SF) 
The SF describes a distribution of momenta and removal 
energies of nucleons inside nucleus. In theoretical mod¬ 
els, the SF combines contributions from the shell model 
(the mean field part) and from the short-range correla¬ 
tions. The second one (about 20% of the overall strength) 
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accounts for the higher momentum and removal energy 
values @,0. Recentlyjthis part has been directly ob¬ 
served experimentally Q. Computations based on the 
SF reproduce precisely the electron scattering data in 
the region of the quasielastic peak 0 . 

In this paper we adopt the plane wave impulse approx¬ 
imation (PWIA) and neglect the effects of final state in¬ 
teractions (FSI) except those from Pauli blocking. Eval¬ 
uation of FSI is necessary in order to compare the predic¬ 
tions of the model with experimental data. It is argued 
in Refs. that in computation of the inclusive cross 

section (after summing over the final hadronic states) the 
construction of the SF and the treatment of FSI are two 
independent problems. Calculation of FSI effects is usu¬ 
ally based on the complex optical potential. In the case 
of electron scattering, FSI do not change much the over¬ 
all cross section at the fixed angle: the elastic peak is 
quenched and the cross section is redistributed making 
the peak wider. For the GeV neutrino interactions, i.e. 
when most of the quasielastic cross section corresponds 
to small values of it is useful to distinguish Pauli 
blocking from the remaining FSI effects, because impact 
of the latter on the cross section is quite small (compare, 
for example fig. 14 in i). 

There is a lot of discussion in the literature on the 
FSI effects in other approaches to quasielastic neutrino- 
nucleus scattering. In some of the recent papers (e.g.. 
Ref. 0), FSI is implied to be responsible only for a re¬ 
distribution of the inclusive cross section among various 
exclusive observable channels. A similar meaning of FSI 
is adopted in many MG generators of events [If . However 
in other theoretical approaches, the FSI effects modify 
also the inclusive cross section. In [l^, it is achieved 
by using the dressed nucleon propagator with the contri¬ 
bution from the self-energy. References [l^ and [I^ de¬ 
scribe the final nucleon wave function as a solution of the 
Dirac equation with the relativistic optical potential, the 
inclusive cross section is then obtained if the imaginary 
part of the potential is neglected. For 1 GeV neutrinos 
the FSI effects are reported to reduce the inclusive cross 
section by ~ 10%. 
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We confine our consideration to quasielastic reactions 
(the primary vertex), but we would like to stress that its 
extension to the resonance region is straightforward. Sev¬ 
eral papers consider the two dynamics togethe r fO . in, 
also in the formalism of the spectral function (Re¬ 
cent discussion of the resonance region can be found 
in [I3,[Il,n') However, we focus on the construction of 
the spectral function, and the inclusion of the resonance 
region is not likely to provide an extra input. In the 
discussion of the resonance region some new theoretical 
uncertainties appear: the properties of resonances in the 
nuclear matter and the nonresonant background [^ . It 
is well known that the theoretical understanding of the 
basic dynamics is still not complete since in the case of 
electron scattering it is difficult to reproduce the cross 
section in the dip region between the elastic and A reso¬ 
nance peaks [2l[. 


For a few-GeV neutrino one has to consider also more 
inelastic channels. They are described by means of the 
DIS formalism with the specific techniques to include 
the nuclear effects. A recent comprehensive study dis¬ 
cussed the following nuclear effects: Fermi motion and 
nuclear binding (FMB), off-shell corrections (OS), nu¬ 
clear pion excess and coherent processes [l^. The re¬ 
sults of the analysis show that in the kinematical region 
of large Bjorken x variable only FMB and OS effects are 
relevant. Both FMB and OS corrections are present in 
our approach: FMB is the core of the spectral function 
computations and OS were included in the de Forest pre¬ 
scription of evaluating weak current hadronic matrix el¬ 
ements [^ . 


There has been a lot of effort to calculate the spec¬ 
tral functions Bam. Systematic computations exist 
only for light nuclei up to oxygen and also for infinite 
nuclear matter. The aim of this paper is to calculate 
the spectral function for argon and to discuss its applica¬ 
tion to computation of the neutrino-argon cross section. 
Liquid-argon technology has been successively tested in 
the case of the ICARUS T600 module and is now 
considered as an option in many neutrino experiments 
[ 2 ^. One of the experiments is the T2K where liquid-Ar 
detector is planned to be an element of the 2 km detector 
complex 271. 

In our computations, we follow the approach proposed 
by Kulagin and Petti which is partially based on the 
paper of Ciofi degli Atti and Simula Q. The computa¬ 
tions contain several approximations, and it is necessary 
to understand how precise the results are. For this reason 
we first analyze the approximation in the case of oxygen, 
the heaviest nucleus for which detailed computations of 
the spectral function exist. Our conclusion is that the 
agreement is satisfactory, and we pursue the computa¬ 
tions for the argon. 


Our paper is organized as follows. In Sec. Ill Al we intro¬ 
duce basic definitions and formulas. Sec. Ill Bl shows how 
the spectral function corresponds to the Fermi gas model 
and in Sec. Ill Cl we present the approach of Kulagin and 
Petti [ 23 . In Sec. [mAI we discuss in detail the pre¬ 


dictions of the spectral function constructed for oxygen 
and compare them with the results obtained by Benhar 
et al. B). Most important are results contained in Secs. 
IIIIBltlnC] In Sec. IIIIBI we construct the argon spec¬ 
tral function. Then, having in mind implementations in 
MC event generators, in Sec. IIII Cl we investigate a pro¬ 
cedure to extract an essential information contained in 
the spectral function. The information is encoded in two 
functions. The first one is the distribution of the nucleon 
momenta, and the second is the average binding energy 
for a given value of momentum. These two functions can 
be inserted in the Smith-Moniz formula for the inclusive 
cross section which can be easily implemented in a MC 
routine. Sec. m contains the conclusions. 


II. SPECTRAL FUNCTION IN PLANE WAVE 
IMPULSE APPROXIMATION 

A. Neutrino-nucleus inclusive cross section 

We denote the initial (final) lepton four-momentum as 
k (k') and the hadronic initial (final) four-momentum as 
p {p^). In general, several hadrons or nuclei fragments 
can be in the final state \f{p^))- The four-momentum 
transfer is q = {to, q) = k — k'. 

The neutrino-nucleus cross section is calculated in 
the Fermi theory approximation (in the neutrino energy 
range of a few GeV the condition is satisfied). 

The inclusive neutrino-nucleus differential cross section 
reads: 


da 


G'jp cos^ 9c 




( 1 ) 


where the leptonic tensor is defined as 

= iTr[7^(l-75)(|f-f TOi)7;.(l-75)(|f'-f mj)] 

= 2{k^kl + k'^kv - k ■ k'- ie^vpak^k"^). (2) 
and the nuclear tensor is 


^ Y,{f{p^)\j>^{o)mA)){f{p^)\j-{o)mA)r 

x{2nfS^{q+PA-p^), (3) 


It is convenient to express the cross section in terms of 
energy and momentum transfer: 


da G% cos^ 6>c |q| ^ 
du!d\q\ 47r E'^ 


( 4 ) 


In this paper we adopt PWIA, i.e. the approximation 
in which neutrino interacts with separate nucleons and 
the nucleon produced in the interaction leaves the nucleus 
without further reinteractions. We assume that 


l/(/)) = \RiPR)) ® \p')- 
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The spectral function -P(p, E) is defined as the proba¬ 
bility distribution to remove from the nucleus a nucleon 
with momentum p and leaving the residual nucleus with 
energy Er = Ma — M + E: 

P{p,E) = Y,^iMA-ER-M + E) 

R, i 

x\{R{pR)\a{p)\i{MA))\'^. (5) 

In our notation, \i{MA)) is the state of the initial nu¬ 
cleus (assumed to be at rest) of mass Ma, \R{pr)) is 
the final nucleus state (after a nucleon of momentum p' 
was ejected) of four-momentum pR = (Er,pr), and M 
denotes the nucleon mass. Summation over all the final 
states and averaging over the spin states of the initial 
nucleus is performed. We omitted spin and isospin de¬ 
pendencies of the annihilation operator, but one has to 
keep in mind that there are two separate spectral func¬ 
tions: one for protons and another for neutrons. 

The spectral function contains a lot of information 
about the nucleus. It follows from its definition that 

nip) = J P{P,E)dE = ^\{R{pR)\a{p)\i{MA))\'^ 

H, i 

= Xl^**^^'4)|a^(p)a(p)|*(M^)) (6) 

i 

is the distribution of momenta inside the nucleus. 

We normalize the nucleon distribution of momenta 
n(p) (and then also the spectral function) as 


adopt the standard de Forest prescription [2^: to ap¬ 
proximate off-shell kinematics one can use free spinors 
and free form factors, taking into account that only a part 
of the energy is transferred to the interacting nucleon, 
and the rest is absorbed by the spectator system. It 
means that one makes a substitution 


q= {uj,q) ^ q = {uj,q) 

(11) 

where u = co — e_B, so that 



(12) 


The model is finally specified by the form of the 
hadronic tensor: 




M2 - -9 H, + ^H,+^e 


'1^ 


H 4 -j- 




.qtJ-p^ 


2M2 


i?5 


(13) 


where 


Hi = Fl{l + t) + t{Fi + F2 f, 

H2 = FI+F^ + tF^, 

H 3 = 2 FAiFi+F 2 ), 

Hi = ^Fiil-T) + ^FiF2 + FAFp-TFl. 

H 5 = H 2 , (14) 

with the notation r = —q‘^j{AM‘^). 


d^pn{p)= / d^pdEP{p,E) = 


Z for protons, 

N for neutrons. 

(7) 

In PWIA it is straightforward to obtain the following 
expression for the inclusive cross section: 


B. Fermi Gas model 

In the FG model 

P{p,E) oc 0(j)F - IpI) 


(15) 


do- ^ G%cos^9c |q| f ^3 -P(p, FI) 
du;d\q\ Att El J ^ E^E^, 

x5{uj + M-E-Ep,)L^,,Hi^r ( 8 ) 

where 




(27r)6 ^ FipAp, (pV'^(O)Ip) (p'IJ'^(O)Ip)* 

spins 



V 2M 


7or‘''''7o 


f + M 
2M 


) 


(9) 


and one only has to decide about the way in which the 
energy is balanced. 

In the previous section we used the energy conservation 
expressed as [compare the 6 function in Eq. ([3])] : 

Ma -\- uj = Epi \J(Ma-\)^ + . (16) 

We recall that the initial nucleus state is assumed to be 
factorized into the neutron with momentum p and the 
spectator system with momentum —p. The neutron en¬ 
ergy can be written as 


and 

= 7'^[7"i(9") + 7 ^ 2 ( 9 ^)] 

+7^5FA(g^)+75^-^^. (10) 

The hadronic tensor H^'^ has to be further specified 
because nucleons described by m are off shell. We will 


\/M2 -I- p2 - esip) = Ep- esip) 

with binding energy eB(p) > 0, so that 

Ma = V{Ma-i)^ + p2 + Ep - eB(p) (17) 
and the energy conservation takes form 

Ep — cb (p) + lo = Epi. 


(18) 
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The Smith-Moniz approach to the FG model is 
to approximate e_B(p) by the constant average value es, 
that is, 


Ep — eB + u! — Epi. 


(19) 


The cross section calculated within this model (with 
Pauli blocking taken into account) is equal to 


d\q\(fpe{pF - IpI) 


das-M Gf 3iV 

dE^ 47rE^ 47rp3 

x6(u; + Ep -€B - Ep') 6»(|p + q| - pf) 

|q' 


X ■ 


EpEpi 




( 20 ) 


Using the concept of the spectral function one can say 
that the Smith-Moniz Q version of the FG model is de¬ 
fined as 



IpI (MeV/c) 

FIG. 1: Momentum distribution of nucleons in [^. 
Dashed and dotted line shows the mean field and correla¬ 
tion contribution, respectively. The solid line represents their 
sum. 


^’s-m(p, U) = M5{\/— CB — M + E) 

x6»(pf- IpI) (21 ) 

where 1b is an average binding energy. 

In the Bodek-Ritchie version of the FG model Q , one 
avoids the approximations in the energy delta function, 
obtaining 

Pb-k{v,E) = N5{Ma - v'(Ma-i) 2 + p2 _ m + U) 
x6»(pf- IpI). (22 ) 

In Ref. 3 , the quadratic momentum distribution is mod¬ 
ified by adding higher momentum tail. A/” is a normal¬ 
ization constant (the same in both cases). 

In ((^ and ((^ Pauli blocking (PB) was ignored, be¬ 
cause the usual PB term present in the expression for 
the cross section [i.e. 0(|p -I- q| — pf)] depends on the 
momentum transfer and therefore is not defined merely 
by characteristics of the nucleus. It is natural to treat it 
separately as a FSI effect. 

C. Spectral function 

In theoretical models, the SF consists of isoscalar and 
isovector parts 

P(p,U) =Psc(p,U) + Pvec(p,i;). (23) 

It is common [^, to think of the isoscalar SF as 

itself consisting of two parts: the mean field (MF) and 
the short-range correlation parts, that is, 

Pscip,E) = [Pmf(P,P) +Pcorr(p,P)] , (24) 

with the normalization given by the condition 

J dE d^p [Pmf(p, E) + Pcorr(p, E)] = 1. (25) 


As a consequence, the momentum distribution ([6]) can 
also be represented as the sum of two contributions: 

»^(p) = —^— Fmf(p) + ncorr(p)J • (26) 

Both nMF(p) and ncorr(p) for a few nuclei are given in [^. 
In Fig. [T] we see how they contribute to the overall mo¬ 
mentum distribution of These data will be used in 
the next section. 

The mean field part Pmf(p,P) provides proper de¬ 
scription of low-energy and low-momentum nucleons: 
they behave as Fermi gas in self-consistent mean field po¬ 
tential with energy levels given by the shell model. The 
presence of the correlation term Pcorr(Pj A) reflects the 
fact that in the ground state A^A^-correlations generate 
high momenta and high removal energies. 

The low-energy contribution can be expressed as [1^ 

Pmf(p,P) = ^ ^ A,n„(p)5(P„ + PF(p)-P) (27) 

OL<i(y.p 

with the recoil energy of the residual nucleus Pf(p) = 
p^/(2M^_i). In above equation A = -|- Z is the total 

number of nucleons, and number of nucleons in the 
shell-model state a of energy Ea and momentum density 
distribution Ua (p) given by momentum space wave func¬ 
tions. In our convention, the values of Ea are positive 
below the Fermi level q:f- 

When we replace Ea by the separation energy E^E av¬ 
eraged over single-particle nucleon levels, the S function 
can be taken outside the sum and the expression simpli¬ 
fies to 

Pmf(p,P) = d{E^^^ + Er{p) - E)^ ^ AaUaip) 

<y.<i(y.p 

= nMY{v)S{E^^'^+E r{p)-E). (28) 

From (e, e'p) experiments we know that (1281) gives, in 
fact, only a part of P(p, E) because the depletion of 
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states is observed of a typical value ^ 0.2. It means 
that ~ 20% of the strength of the spectral function must 
be represented by a physics outside the shell model. It 
is usually explained in terms of short-range correlations 
producing pairs of high-momentum nucleons. To model 
Id:orr(p, E), One assumes that two nucleons form a cluster 
with high relative momentum and other (A — 2) nucle¬ 
ons remain soft [2^. Within this picture the following 
formula can be obtained 12211: 


/ X / ^ ^ /“ 

Ecorrip^ E) = ncorrlP)] 7\ 


X [exp(-apL„) - exp(-apLx)] > (29) 

where a = 3/(4(p^)/3) with the mean square of momen¬ 
tum (p^), P = {A — 2)/{A— 1), and 


P^ ■ = 

Ir'min 


— 


P\p\ - V2Mp[E-E<^^^ -Er{p)] 
P\p\ + ^2Mp[E - - £;fi(p)] 


(30) 


The constant E^^'^ is the two-nucleon separation energy 
averaged over configurations of the (A—2) nucleon system 
with low excitation energy. 

Even in nuclei with equal number of neutrons N and 
protons Z there is a nonzero isovector contribution com¬ 
ing from Coulomb interaction and isospin-dependent ef¬ 
fects in other interactions, but the standard approach is 
to neglect differences between the neutron and proton SF. 
Taking into account differences between the neutron and 
proton energy levels and their occupation probabilities in 
the approach presented here yields only small change of 
the threshold energy E^^\ so it only slightly shifts the 
plot of the differential cross section. The reason of us¬ 
ing the proton momentum distribution in both cases is 
purely practical: it’s known with better accuracy. 

For nonsymmetric nuclei, such as argon |gAr, the neu¬ 
tron SF must be different than the proton one. Due to 
the Pauli exclusion, the surplus neutrons can occupy only 
the Fermi level and thus it is assumed that the dominant 
contribution to Pvec{p,E) comes from momenta and en¬ 
ergies very close to the Fermi level. Kulagin and Petti 
propose the approximation 

Pvec(p, E) = 5(|p| - Pf)S(E - Ef). (31) 

2 A-kpf 


III. RESULTS 


TABLE I: Occupation probabilities Cni and numbers of 
particles Z^i ~ \Ani for proton shell-model states in ^® 0 . 


State 

160 

Znl 

2 s 

0.06 

2 

Id 

0.14 

10 

OiF 



ip 

0.80 

6 

Is 

0.87 

2 


TABLE II: Energy spectrum of protons and neutrons in 
[sJl (the data for lsi /2 and ap are taken from [3^1. 


State a 

Protons 

Neutrons 

ld3/2 

-4.65 

-0.93 

2Si/2 

0.08 

3.27 

ld5/2 

0.59 

4.15 

olf 


10 

lPl/2 

12.11 

15.65 

lP3/2 

18.44 

21.8 

lSl /2 

45 

47 


One-nucleon separation energy E^^'> has been intro¬ 
duced in (l28ll as an average energy of the single-particle 
states below the Fermi level, so 

2^ _^ r 2 _^ 

E^^'> = ^ X! / ’^a(p)d^P=^ ^ AaEaCa, 

Q;<Q!jr Q!<aj7- 

(32) 

where Cq, is the occupation probability of the state a. 

Multi-index a is a compact way of writing quantum 
numbers (nlj); therefore A^ = Anij, etc. The occupation 
probabilities Cnij for oxygen were not available to us, 
but we found values of Cni defined as 

2(1 -f ^) -f 1 ^ 2(1 — ^) -|- 1 

“ 2(21 -b 1) '^"■'('+ 5 ) 2(21 -b 1) 

^ + ^ ^ /oo\ 

“ 2^+1 2r+T 

It follows from the data for calcium in Sec. IIIIBl that the 
values of Cnij for the levels differing only in spin state 
are very similar. Assuming that c„ 1 ( 1 + 1 ) ~ (or 

equivalently that ~ Cni and « Cni) and 

using data from Tables III and im we get the following val¬ 
ues: 


A. Verification of the method - an example of 
oxygen 


EiIq = 19.18 MeV for protons, 
E’fgQ = 22.32 MeV for neutrons. 


It is important to check how good the presented ap¬ 
proximation is for oxygen, for which detailed computa¬ 
tion of the spectral function exists Q. To perform the 
comparison we have to calculate the values of the con¬ 
stants appearing in ((^ and (1^ : E^^\ E^‘^\ and (p^). 


The Benhar’s spectral function was calculated for pro¬ 
tons; therefore, to make appropriate comparisons, we 
used e[Iq = 19.18 MeV in our numerical calculations. 
We checked also that applying the neutron value does 
not change the cross sections significantly. 
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FIG. 2: Comparison of the quasielastic differential cross 
section da/dE^ of for the Benhar’s spectral function 
(dashed), our spectral function (solid), and the Fermi gas 
model with pp = 225 MeV (dotted line). Pauli blocking is 
not included. 




FIG. 4: Comparison of the nucleon momentum distribution 
computed from the Benhar’s spectral function of and the 
distribution given by Ciofi degli Atti and Simula used in 
our computations. 



FIG. 3: Same as Fig. [2] but with Pauli blocking (details in 
the text). 

Two-nucleon separation energy can be approxi¬ 
mated by 


FIG. 5: Contributions to the quasielastic differential cross 
section da/dEf^ of in our spectral function with Pauli 
blocking (jip = 209 MeV). Mean field and correlation con¬ 
tribution are shown by dashed and dotted line, respectively. 
The solid line represents their sum. 


= Ma-2 + 2M-Ma, (34) 


and for oxygen we obtained 

= 26.33 MeV. 

The mean value of is 

^ /p^nMF(p) (Pp 
^ / nMF(p) d?P 


(35) 


Using the mean field momentum distribution for as 
in we obtain 


(p2)ieo = J(fic)20.9 fm"2 = (162.2 MeV/c)^. 

Figure [5] shows a comparison of predictions for the oxy¬ 
gen differential cross section in energy transfer given by 


our spectral function with predictions given by the Ben¬ 
har’s spectral function and the Fermi gas model. All 
numerical results in this paper have been obtained using 
the dipole parameterization of the form factors and axial 
mass Ma = 1.03 GeV. 

To make comparisons with other numerical computa¬ 
tions, we included the Pauli blocking effect, taking the 
average Fermi momentum calculated from the known 
density profile of a nucleus. The values we use are 
Pp = 209 MeV 0 and pp = 217 MeV for oxygen and 
argon, respectively. Note that the corresponding values 
needed for the FG model to reproduce the quasielastic 
peak in electron scattering (namely, pp = 225 MeV 
and Pp = 251 MeV) are visibly higher. We checked that 
the use of the same value for Pauli blocking in both FG 
and SF does not reduce discrepancy between them sig¬ 
nificantly. 

The difference between predictions of two spectral 
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functions is not large, and it is greater if the Pauli block¬ 
ing is included, see Fig [31 We identify two reasons why 
the discrepancy is observed. A part of it is explained 
by the different momentum distributions in the spectral 
functions (see Fig. II. We deduce that the difference 
must be also caused by oversimplified treatment of the 
Fmf(p, E)^ because for high this part gives dominant 
contribution to the differential cross section, as shown in 
Fig. |5l 

In spite of simplicity of our approach, the results are 
in satisfactory agreement with the “exact” spectral func¬ 
tion, and we continue our investigation to obtain a de¬ 
scription of the argon nucleus. 


TABLE III: Energy Ea, occupation probability c^, and num¬ 
ber of particles Za = for proton shell-model states in 

IgCa 


a 

Eol 

Ca. 

Eol 

157/2 

-21.00 

0.04 

8 

159/2 

-9.75 

0.06 

10 

1/5/2 

-4.51 

0.08 

6 

2pi/2 

-2.02 

0.07 

2 

‘ 2 P 3/2 

-0.92 

0.08 

4 

1/7/2 

1.15 

0.14 

8 

qf 

4.71 



1^3/2 

8.88 

0.85 

4 

2Si/2 

10.67 

0.87 

2 

1^5/2 

14.95 

0.88 

6 

lPl/2 

31.62 

0.91 

2 

lP3/2 

36.52 

0.92 

4 

lSl/2 

57.38 

0.93 

2 


TABLE IV: Energy Ea, occupation probability Ca, and num¬ 
ber of particles Na = ^Aa for neutron shell-model states in 
IgGa 0. 

a 

Eol 

Ca 

Na 

1 / 5/2 

1.50 

0.07 

6 

2Pl/2 

4.19 

0.07 

2 

2P3/2 

5.59 

0.08 

4 

1 / 7/2 

8.54 

0.12 

8 

olf 

12 



1^3/2 

15.79 

0.89 

4 

2si/2 

17.53 

0.91 

2 

1^5/2 

22.48 

0.93 

6 

lPl/2 

39.12 

0.96 

2 

14 * 3/2 

43.8 

0.96 

4 

lSl/2 

66.12 

0.97 

2 


B. Argon spectral function 


In nature. 

the most common isotope of argon is 

fOAr. 

The feature of this nucleus which was not present : 

in the 

previous calculations is 

its lack of isospin symmetry. It 

is obvious that for such nuclei, the vector contribution to 



FIG. 6: Quasielastic differential cross section da/dEfj, off igAr 
as a function of produced muon energy for the Fermi gas 
model {pf = 251 MeV, es = 28 MeV) and the spectral func¬ 
tion. Pauli blocking is absent. 



FIG. 7: Same as Fig. (6] but Pauli blocking is present. 


the spectral function cannot be neglected, and therefore 
we estimate it using Eq. m- 

We lack data on the energy levels and occupation prob¬ 
abilities for the argon. The best we can do is to assume 
that they do not differ significantly from those for cal¬ 
cium given in Tables IIIII and IIVI and use them to obtain 
the value of one-nucleon separation energy. In this way 
we get for neutrons 

= 29.26 MeV. 

From Eq. (j34p two-nucleon separation energy for argon 
is 

Eiol = 16.46 MeV. 

We approximate the momentum distribution of argon 
with the one of calcium [^, what yields 

(p')«Ar = lihcfim W = (177.2 MeV)2. 

Using the argon spectral function, we calculated the 
inclusive cross section for neutrino-argon interaction 

























FIG. 8: Contributions to the differential cross section shown 
in Fig. [3 from PMF(p,iF), Pcorr(p,£') and Pvec(p,P). 



FIG. 10: Quasielastic differential cross section dcr/dP^ off 
obtained from the effective description, the Benhar’s spectral 
function and the Fermi gas. Pauli blocking is present. 


first without (Fig. [6]) and then with the Pauli block¬ 
ing (Fig. [7]). In Fig. [51 we show the contributions from 
P^ec{p,E), Pmf{p,E) and Pcorr(p,£^)- We see that 
Pvec(p, P) is clearly responsible for the singular behav¬ 
ior of the differential cross section at small values of the 
energy transfer. The behavior at the quasielastic peak 
is determined mainly by the Pmf(p, P) and the detail 
knowledge of argon energy levels can lead to modifica¬ 
tions in this region. We also compared the differential 
cross section per neutron for oxygen and argon and found 
that they are rather similar. 


C. Effective description of nucleus 


complicate model too much. This approach is simpler 
than the one based on the complete spectral function 
and is easy to apply in MC simulations. We refer to it as 
the effective description. It is defined by two functions 
eB(p) and ?T,(p) which characterize specific nucleus and 
can be calculated directly from the spectral function 

nip) = J dEP{p,E) 

and 

esip) = ^ J dE{^/M^ + p^-M+E)Pip,E). (36) 


Many Monte Carlo event generators rely on the Smith- 
Moniz version of the FG model despite the fact that 
the energy conservation given by Eq. (HU) and the step- 
function momentum distribution provide a rather crude 
approximation of the real nucleus. 

The use of the energy conservation written in ITSl) and 
the momentum distribution of a specific nucleus does not 



IpI (MeV/c) 

FIG. 9: Binding energy dependence on nucleon momentum 
^^(p) extracted from the Benhar’s spectral function. 


The momentum distribution n(p) obtained from the Ben¬ 
har’s spectral function for oxygen was presented in Fig.jU 
and the binding energy dependence on momentum is 
shown in Fig. jS] One can see that high-momentum nu¬ 
cleons are deeply bound [2^ and they cannot leave a 
nucleus. One can also see that for |p| £ [0;200] MeV, 
the binding energy is roughly constant, varying between 
~ 35 and ~ 47 MeV. 

The cross section in the effective description is of the 
form 


dfJefF 

dE^ 


Gp cos^ dc 
X 5 ioj Ep 


N f 

<iiqi<iy»(p) 

-eB(p) - Ep/) 6»(|p-hq| -pp) 


X 


|q| 

EpEpi 


Lp^E[)^g. 


(37) 


The only difference between Lp^El^^ and lies 

in the energy transfer, because now uj = to — eB(p). 

Figs. [TUHm show results obtained within this model. 

The effective description of oxygen makes use of the 
Benhar’s spectral function. The cross section daes/dEp 
(Fig. [TUll is only slightly different with respect to the one 
calculated with the spectral function. The difference is 
seen at low energy transfer. We notice that the effective 
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FIG. 11: Quasielastic differential cross section dcr/d|p'| off 
obtained from the effective description, the Benhar’s 
spectral function and the Fermi gas {pf = 225 MeV, eB=27 
MeV). Pauli blocking is absent. 



FIG. 12: Quasielastic differential cross section in = —q^ 
of obtained from the effective description, the Benhar’s 
spectral function and the Fermi gas {pf = 225 MeV, Is=27 
MeV). Pauli blocking is present. 


description represents a major improvement over the FG 
model. 

We compare also predictions of the three models for the 
distribution of nucleons ejected from oxygen. The results 
are shown in Fig. 1111 In the computations we did not 
apply the Pauli blocking. Its straightforward implemen¬ 
tation simply eliminates nucleons with momenta lower 
then the Fermi momentum. Predictions of the spectral 
function and the effective description are almost identi¬ 
cal. It is because the distributions of momenta of two 
approaches are the same. We see again that the effec¬ 
tive description is a good approximation of the spectral 
function. 

As is shown in Fig. na the effective model also succeeds 
in precisely reproducing the SF’s differential cross section 
daldQ"^. 

Finally in Figs. [T3] and [H] we show a comparison of 
predictions from our spectral function of argon with the 
effective description based on it. The differences between 



FIG. 13: Quasielastic differential cross section da/dE/j, of ig Ar 
obtained from the effective description, our spectral function 
and the Fermi gas {pf = 251 MeV, es = 28 MeV). Pauli 
blocking is present. 



FIG. 14: Quasielastic differential cross section d(T/d|p'| of 
igAr obtained from the effective description, our spectral 
function and the Fermi gas {pf = 251 MeV, Is = 28 MeV). 
Pauli blocking is absent. 


the two models are smaller then in the case of oxygen 
because of the way in which the dominant mean field 
part of the argon SF was treated. We notice that the 
singularity present in Fig. [13] caused by the treatment 
of isovector part of the spectral function is not seen in 
Fig [Til 


IV. CONCLUSIONS 

We constructed the spectral function for argon and ap¬ 
plied it in the PWIA computation of the CC quasielastic 
cross section of neutrino-argon scattering. The construc¬ 
tion method was based on the ideas contained in [2^ . We 
verified that when applied to oxygen, the method leads 
to the results similar to those obtained from the more 
detailed spectral function obtained by Benhar [^. The 
only FSI effect included in the discussion was Pauli block¬ 
ing. We proposed also the effective approach in which the 
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information contained in the spectral function is used to 
define two functions. We argued that this approach yields 
results similar to the spectral function, but at the same 
time it can be easily applied in MC routines. 

We see two possible applications of our results. The 
argon spectral function can be used in more precise com¬ 
putations of the neutrino cross sections in a few GeV 
energy region. The effective approach may be used in ex¬ 
isting MC event generators giving rise to the important 
improvement with respect to the Fermi gas model. For 
example, Fig. [T^] shows that it reconstructs well the 
distribution of events which is important because of the 
observed low deficit of events in the K2K experiment. 


The effective approach can be also used in cross section 
computations for other targets important in neutrino ex¬ 
periments like carbon and iron for which appropriate SFs 
exist 
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